Frequency Dependence of Magnetopolarizability of Mesoscopic Grains 
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We calculate average magnetopolarizability of an isolated 
metallic sample at frequency ui comparable to the mean level 
spacing A. The frequency dependence of the magnetopolar- 
izability is described by a universal function of w/A. 

PACS numbers:73.20.Fz,73.23.-b,73.61.-r 



Mesoscopic effects in electric polarizability of small 
metallic particles have been discussed in the literature 
starting from the seminal work by Gor'kov and Eliash- 
berg (GE) (l). The role of screening effects, which were 
not taken into account in the original paper Q , has been 
subsequently emphasized |Q. Whereas GE predicted a 
giant effect of level correlations on the polarizability, 
screening actually reduces this quantum effect to a rela- 
tively small correction to the classical value. Like other 
mesoscopic effects related to quantum interference, this 
correction is affected by magnetic field. It can thus be 
observed in magnetopolarizability of a grain, similarly 
to the weak localization correction to resistivity 
An experimental study of quantum corrections in elec- 
tromagnetic response of an ensemble of 10 5 mesoscopic 
samples has been reported recently ||[JJ . 

The aim of this paper is to calculate the dependence 
of the magnetopolarizability on the frequency to of the 
applied electric field in the range ui ~ A, where A is 
the mean level spacing. This dependence is given by a 
universal function of the parameter w/A. Note that the 
frequency in the experiments [[H?]] was comparable to A, 
and thus experimental verification of this scaling form is 
in principle feasible. 

Previously we have demonstrated [^|J5| how the su- 
persymmetric sigma-model can be employed to calculate 
magnetopolarizability in the grand canonical ensemble. 
However, the experimentally relevant situation of isolated 
grains is appropriately described in the framework of the 
canonical ensemble. The crucial role of the statistical en- 
semble for the problem considered has been emphasized 
in Ref. |J (see also Ref. Technically, calculation 

of the dynamical response within the canonical ensemble 
is more complicated, since it involves three-level corre- 
lation functions, which poses a serious obstacle to the 
application of the supersymmetry technique. Below, we 
bypass this difficulty by combining information on statis- 
tical properties of wavefunctions, obtained with the su- 
persymmetry technique, with the three-level correlation 
function known from the random matrix theory (RMT). 



The (complex-valued) polarizability is defined as 
d(uS) — a(u>)E(u>), where E and d are the external elec- 
tric field and the induced dipolc moment, respectively. 
The expression for the ensemble-averaged magnetopolar- 
izability = a(u>,B) — a(u>,0), where B is an ap- 
plied magnetic field (strong enough to break the time- 
reversal symmetry) has the form |],|5| 
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where <£>(r) is the local potential in the sample (resulting 
from screening of the external electric field), the symbol S 
denotes the difference between quantities with and with- 
out magnetic field, and II is the polarization operator, 

II(ri,r 2 ;w) = ( $Z ^( r i)^n(riW n (r 2 )^ m (r 2 ) 



(2) 



Here m and n label exact single-particle states, the an- 
gular brackets denote the ensemble averaging, and it is 
assumed that u> has an (infinitesimally small) positive 
imaginary part, u> = u) + iO. We will consider the low- 
temperature limit, T <C A, thus setting T = in the 
sequel. Substitution of (0) into (Tr!) yields 
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where $ m „ = / dril)* m {r)<5>(r)yb n {r). 

The canonical ensemble is realized by pinning the 
Fermi-level to one of the single-particle levels e F — 
f-k + 0. Splitting the sum in (|^) into two contributions 
with n = k and n ^ k, we find 
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where i?2(e) and ^(e, ei) are the two-level and the three- 
level correlation functions (normalized to unity at e, e\ 3> 
A), 
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R 2 (e) = A 2 (^S{E-E i )S(E + e-E j )J ; 
R 3 (e, ei) = A 3 (Y / 5(E- E t )S(E + e - E 3 ) 

\ ijk 

x 5{E + Cl - E^ , (5) 
and (| < i ) | 2 ) £ is the average squared matrix element, 
(|$| 2 > e = J R 2 " 1 (e)A 2 /^|$ i ,| 2 



x 5(E - Ei)S{E + e-Ej 



(6) 



When writing Eq. (Q), we decoupled the wave function 
correlations (^) from the three-level correlation function 
i?3(e, ei). Indeed, we know from the supersymmetry cal- 
culations of two-level correlation functions |ll|] that (i) 
the level correlation function has the RMT form, up 
to l/g 2 corrections; (ii) the wavefunction correlations 
have no dependence on lo/A in the order l/g. Here 
g ~ E c /A ^> 1 is the dimensionless conductance of the 
grain, E c is the Thouless energy, and we consider the 
frequency range cj <C E c . We make thus an (extremely 
plausible) assumption that these properties hold also for 
higher order correlation functions, which allows us to pro- 
ceed in the case of the canonical ensemble. 

We begin the evaluation of Eq. ([|) by considering the 
term which contains the three-level correlator. We write 
i? 3 (e, ei) = R 2 (e) + Ra{e, ei) and denote the correspond- 
ing contributions as and . In the leading order in 
l/g the e-integral in is a sum of contributions from 
the regions e ~ E c and e ~ A (to be denoted as ct% a ^ 

and respectively). In the former we can neglect ui, 
which yields after taking into account the_orthogonality 
and completeness of the eigenfunctions 
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where we defined 
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<l*| 2 )o=(l (I) a« £c (8) 
= 72 / ^i^2*(T"i)*(r2)nr,(ri,r 2 ) , 

and Il£> is the diffusion propagator satisfying 

- DV 2 n D (n,r 2 ) = {ttu)- 1 [8(n - r a ) - V' 1 ] , (9) 

with the boundary condition V^II^ = 0. In the con- 
tribution otg we can neglect the e-dependence of the 
matrix element, thus replacing it by (|<I>| 2 ) , 
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The same is valid for the term a B , as well as for the re- 
maining contribution of the first term in square brackets 
in (0). Collecting everything, we finally get 
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where 



F(lo) = 1 + 
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Since the integrals in (|T^) are determined by the range 
e, ei ~ A, we can use the RMT results |12 for the level 
correlation functions R 2 and R3 entering this formula. 
Therefore, F(u>) is in fact a universal function of the di- 
mensionless parameter s = lo/A. Real and imaginary 
parts of this function describe the influence of the mag- 
netic field on polarizability and absorption of grains, re- 
spectively, and are plotted in Fig. 1. 




FIG. 1. Real and imaginary parts of the function F [Eq. 
(0)] representing the frequency dependence of the magne- 
topolarizability as (a;). 

Note that F(s) = at s = 0, which is an identity relat- 
ing the two-level and the three-level correlation functions 
entering Eq. ([l2]). This identity can be derived by using 
the invariance of the level correlations with respect to a 
perturbation fl3[| . 

As is seen from Fig. 1, our results predict a positive 
magnetopolarizability for all values of u> and a negative 
magnetoabsorption for almost all lo. This is in qualita- 
tive agreement with the findings of Ref. M, where the 
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measurements have been performed at a frequency cor- 
responding to s ~ 0.21. It would be very interesting to 
have the experimental data for several frequencies in the 
range lo ~ A in order to check our prediction for the 
universal scaling function F (s) describing the frequency 
dependence of the magnetopolarizability induced by the 
level and eigenfunction statistics in the grains. 
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